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We study the synchronization of coupled maps on a variety of networks including regular one and 
two dimensional networks, scale free networks, small world networks, tree networks, and random 
networks. The dynamics is governed by a local nonlinear map for each node of the network and 
interactions connecting different nodes via the links of the network. For small coupling strengths 
nodes show turbulent behavior but form phase synchronized clusters as coupling increases. We 
identify two different ways of cluster formation, self-organized clusters which have mostly intra- 
cluster couplings and driven clusters which have mostly inter-cluster couplings. The synchronized 
clusters may be of dominant self-organized type, dominant driven type or mixed type depending 
on the type of network and the parameters of the dynamics. We also observe ideal clusters of both 
self-organized and driven type. There are some nodes of the floating type that show intermittent 
behaviour between getting attached to some clusters and evolving independently. The residence 
times of a floating node in a synchronized cluster show an exponential distribution. We define 
different states of the coupled dynamics by considering the number and type of synchronized clusters. 
For the local dynamics governed by the logistic map we study the phase diagram in the plane of 
the coupling constant (e) and the logistic map parameter (p). For large coupling strengths and 
nonlinear coupling we find that the scale free networks and the Caley tree networks lead to better 
cluster formation than the other types of networks with the same average connectivity. For most 
of our study we use the number of connections of the order of the number of nodes which allows 
us to distinguish between the two mechanisms of cluster formation. As the number of connections 
increases the number of nodes forming clusters and the size of the clusters in general increase. 
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I. INTRODUCTION 

Several complex systems have underlying structures 
that are described by networks or graphs and the study 
of such networks is emerging as one of the fastest grow- 
ing subject in the physics world One significant 
discovery in the field of complex networks is the obser- 
vation that a number of naturally occurring large and 
complex networks come under some universal classes and 
they can be simulated with simple mathematical models, 
viz small- world networks |3|, scale- free networks [j] etc. 
These models are based on simple physical considera- 
tions and have attracted a lot of attention from physics 
community as they give simple algorithms to generate 
graphs which resemble several actual networks found in 
many diverse systems such as the nervous systems 
social groups [a, world wide web , metabolic networks 
01, food webs 13 and citation networks p^ . 

Several networks in the real world consist of dynamical 
elements interacting with each other. These networks 
have a large number of degrees of freedom. In order to 
understand the behaviour of these systems we study the 
synchronization and cluster formation of these dynamical 
elements evolving on different networks and connected 
via the links of the networks. 
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Synchronization and cluster formation lead to rich 
spatio-temporal patterns when opposing tendencies com- 
pete; the nonlinear dynamics of the maps which in the 
chaotic regime tends to separate the orbits of different el- 
ements, and the couplings that tend to synchronize them. 
There are several studies on coupled maps/oscillators 
on regular lattices as well as globally coupled networks. 
Coupled map lattices with nearest neighbor or short 
range interactions show interesting spatio-temporal pat- 
terns, and intermittent behavior JA, ,12]. Globally cou- 
pled maps (GCM) where each node is connected with 
all other nodes, show interesting synchronized behavior 
[T3| . Formation of clusters or coherent behaviour and 
then loss of coherence are described analytically as well 
as numeri cally at different places with different points of 
views [3|ll[ll|ll|ll[ll. Chaotic coupled map lat- 
tices show beautiful phase ordering of nodes "2^. There 
are also some studies on cou pled maps on different types 
of networks. Refs. 0, 123. l23|| shed some light on the 
collective behavior of coupled maps/oscillators with lo- 
cal and non-local connections. Random networks with 
large number of connections also show s ync hronized be- 
havior for large coupling strengths l25l |2^ . There 
are some studies on synchronization of coupled map s on 
Cayley tree ^2?! , small- world networks |28l |2^ l30 f| and 
hierarchal organization (3l| . Coupled map lattice with 
sine-circle map gives synchronization plateaus [s^ • Ana- 
lytical stability condition for synchronization of coupled 
maps for different types of linear and non- linear couplings 
are also discussed in several papers 0[ |3 113 - Syn- 
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chronization and partial synchronization of two coupled 
logistic maps are discussed at length in Ref. (sj- Apart 
from this there are other studies that explore different 
properties of coupled maps [sS Isglliollillli^ . 

Coupled maps have been found to be useful in several 
practical situations. These include fluid dynamics [43l |. 
nonstatistical behavior in optical systems 44] , convection 
EH l4^- stock markets , ecological systems [i^ , logic 
gates [4^, solitons l50| and c-elegans jfilj . 

Here we study the detailed dynamics of coupled maps 
on different networks and investigate the mechanism of 
clustering and synchronization properties of such dynam- 
ically evolving networks. We explore the evolution of in- 
dividual nodes with time and study the role of different 
connections in forming the clusters of synchronized nodes 
in such coupled map networks (CMNs). 

Most of the earlier studies of synchronized cluster for- 
mation have focused on networks with large number of 
connections (~ N"^). In this paper, we consider networks 
with number of connections of the order of N. This small 
number of connections allows us to study the mechanism 
of synchronized cluster formation and the role that differ- 
ent connections play in synchronizing different nodes. We 
identify two phenomena, driven and self-organized phase 
synchronization js^ . The connections or couplings in 
the self-organized phase synchronized clusters are mostly 
of the intra-cluster type while those in the driven phasen 
synchronized clusters are mostly of the inter-cluster type. 
As the number of connections increases more and more 
nodes are involved in cluster formation and also the cou- 
pling strength region where clusters are formed increases 
in size. For large number of connections, typically of 
the order of and for large coupling strengths, mostly 
one phase synchronized cluster spanning all the nodes is 
observed. 

Depending on the number and type of clusters we de- 
fine different states of synchronized behaviour. For the 
local dynamics governed by the logistic map, we study 
the phase diagram in the /i — e plane, i.e. the plane de- 
fined by the logisting map parameter and the coupling 
constant. 

The paper is organized as follows. In section II, we give 
the model for our coupled map networks. We also define 
phase synchronization and synchronized clusters as well 
as discuss the mechanisms of cluster formation. In sec- 
tion III, we present our numerical results for synchroniza- 
tion in different networks and illustrate the mechanism 
of cluster formation. This section includes the study of 
the phase diagram, lyapunov exponent plots, behavior of 
individual nodes, dependence on number of connections, 
and behaviour for different types of networks. Some uni- 
versal features of synchronized cluster formation are dis- 
cussed in section IV. Section V considers circle map. Sec- 
tion VI concludes the paper. 



II. COUPLED MAPS AND SYNCHRONIZED 
CLUSTERS 

A. Model of a Coupled Map Network (CMN) 

Consider a network of N nodes and Nc connections (or 
couplings) between the nodes. Let each node of the net- 
work be assigned a dynamical variable ,i = 1,2, . . . , N . 
The evolution of the dynamical variables can be written 
as 

N 

xU^ = il~e)f{xl} + ^Y.^,,g{xl), (1) 

J = l 

where xl is the dynamical variable of the i-th node at 
the t-th time step and e is the coupling strength. The 
topology of the network is introduced through the adja- 
cency matrix C with elements Cy taking values 1 or 
depending upon whether i and j are connected or not. 
C is a symmetric matrix with diagonal elements zero. 
ki — ^ij is t-^is degree of node i. The factors (1 — e) in 
the first term and ki in the second term are introduced 
for normalization. The function f{x) defines the local 
nonlinear map and the function g{x) defines the nature 
of coupling between the nodes. Here we present detailed 
results for the logistic map, 

f{x)^fix{l~x) (2) 

governing the local dynamics. We have also considered 
some other maps for local dynamics. We have studied 
different types of linear and non- linear coupling functions 
and here discuss the results in detail for the following two 
types of coupling functions. 

g{x) = X (3) 
9{x) = /(^)- (4) 

We refer to the first type of coupling function as linear 
and the later as nonlinear. 

B. Phase synchronization and synchronized 
clusters 

Synchronization of coupled dynamical systems [ssL Is^ 
IsSi] is indicated by the appearance of some relations be- 
tween the functionals of different dynamical variables due 
to the interactions. The exact synchronization corre- 
sponds to the situation where the dynamical variables for 
different nodes have identical values. The phase synchro- 
nization corresponds the situation where the dynamical 
variables for different nodes have some definite relation 
between the phases [s^ [s^, IH, ■ When the number 
of connections in the network is small ('^ A^) and when 
the local dynamics of the nodes (i.e. function f{x)) is 
in the chaotic zone, only few clusters with small number 
of nodes show exact synchronization. However, clusters 
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with larger number of nodes are obtained when we study 
phase synchronization. For our study we define the phase 
synchronization as follows 

Let Vi and Uj denote the number of times the dynam- 
ical variables x\ and x^, t = t^.to -\- 1,2, . . . ,to + T — 1, 
for the nodes i and j show local minima during the time 
interval T starting from some time t^. Here the local 
minimum of x\ at time t is defined by the conditions 
x\ < xl_i and < xj^j^. Let Vij denote the number 
of times these local minima match with each other, i.e. 
occur at the same time. We define the phase distance, 
da, between the nodes i and j by the following relation 



(5) 



Clearly, u,,^ 



dji = d 



Also, 



dij — 



when all minima of 



variables and x^ match with each other and d, 
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when none of the minima match. In Appendix A, we 
show that the above definition of phase distance satisfies 
metric properties. We say that nodes i and j are phase 
synchronized if dij = 0, and a cluster of nodes is phase 
synchronized if all the pairs of nodes belonging to that 
cluster are phase synchronized. 



C. Mechanism of cluster formation 

Now we consider the relation between the synchronized 
clusters which are formed by the dynamical evolution of 
nodes and the coupling between the nodes of network 
which is a static property for a given network represented 
by the adjacency matrix. Clustering is obviously because 
of the coupling between the nodes of the network and may 
be achieved in two different ways |52j |. 

(i) The nodes of a cluster can be synchronized because 
of intra-cluster couplings. We refer to this as the self- 
organized synchronization. 

(ii) Alternately, the nodes of a cluster can be synchro- 
nized because of inter-cluster couplings. Here nodes of 
one cluster are driven by those of the others. We refer to 
this as the driven synchronization. 

We are able to identify ideal clusters of both the types, 
as well as clusters of the mixed type where both ways 
of synchronization contribute to cluster formation. We 
will discuss several examples to illustrate both types of 
clusters. 



(b) Partially ordered state (HI): Nodes form a few clus- 
ters with some nodes not attached to any clusters. 

(c) Ordered state (IV): Nodes form two or more clusters 
with no isolated nodes at all. This ordered state can be 
further divided into 3 substates based on the nature of 
nodes belonging to a cluster. These partitions are chaotic 
ordered state, quasi-periodic ordered state, and periodic 
ordered state. 

(d) Coherent state (V) : Nodes form a single synchronized 
cluster. 

(e) Variable state (II): Nodes form different states, par- 
tially ordered or ordered states depending on initial con- 
ditions. 

In addition to the above definition of states depending 
on the number of clusters, we further divide the states 
having synchronized clusters into subcategories depend- 
ing on the type of clusters i.e. self-organized (S), driven 
(D) or mixed type (M). 



III. NUMERICAL RESULTS 

Now we present the numerical results of the coupled 
dynamics of variables associated with nodes on different 
types of networks. Starting from random initial condi- 
tions the dynamics of Eq. after an initial transient, 
leads to interesting phase synchronized behavior. The 
adjacency matrix C depends on the type of network and 



1 if the corresponding nodes in the network are 



connected and zero otherwise. First we will discuss our 
numerical results in detail for scale- free network and then 
we will discuss other networks. 



A. Coupled maps on scale-free network 

1. Generation of Network 

The scale free network of N nodes is generated by using 
the model of Barabasi et.al. 63] . Starting with a small 
number, mg, of nodes, at each time step a new node with 
TO < Too connections is added. The probability 'K{ki) that 
a connection starting from this new node is connected to 
a node i depends on the degree ki of node i (preferential 
attachment) and is given by 



ih + 1) 



D. States of synchronized dynamics 

We find that as network evolves, it splits into several 
phase synchronized clusters. The states of the coupled 
evolving system on the basis of number of clusters can 
be classified as in Ref. 62] 

(a) Turbulent state (I) : All nodes behave chaotically with 
no cluster formation. 



After T time steps the model leads to a network with 
N = T + toq nodes and tot connections. This model 
leads to a scale free network, i.e. the probability P{k) 
that a node has a degree k decays as a power law, 

P{k) - /fc^, 

where A is a constant and for the type of probability 
law 7r(fc) that we have used A = 3. Other forms for the 
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FIG. 1: Phase diagram showing different regions in the two 
parameter space of fi and e for scale free network for f{x) — 
fj,x{l — x) and g{x) = x. Different regions based on number 
of clusters are I. Turbulent region, II. region with varying be- 
haviour, III. Partially ordered region, IV. Ordered region, V. 
Coherent region. The symbols T, S, M, D, P, Q and F respec- 
tively correspond to turbulent behaviour, self-organized syn- 
chronization, mixed synchronization, driven synchronization, 
periodic, quasiperiodic and fixed behaviour. Region bound- 
aries are determined based on the asymptotic behaviour us- 
ing several initial conditions, number of clusters and isolated 
nodes, synchronization behaviour and also the behaviour of 
the largest Lyapunov exponent. The dashed lines indicate un- 
certainties in determining the boundaries. Calculations are 
for A'^ = 50, m — 1,T = 100. The inset shows the phase 
diagram for the entire range of parameter fj, i.e. from to 4. 

probability 7r(fc) are possible which give different values 
of A. However, the results reported here do not depend 
on the exact form of 7r(fc) except that it should lead to a 
scale-free network. 



2. Linear coupling 

Phase diagram: First we start with the linear coupling, 
g{x) — X. Fig. ^shows the phase diagram in the two pa- 
rameter space defined by fi and e for the scale-free net- 
work with m — Too = l,N ~ 50, T = 100. For fj, < 3, we 
get a stable coherent region (region V-F) with all nodes 
having the fixed point value. To understand the remain- 
ing phase diagram, consider the line /i = 4. Fig. [21 shows 
the largest Lyapunov exponent A as a function of the cou- 
pling strength e for fi — A. We can identify four different 
regions as e increases from to 1; namely the turbulent 
region, the variable region (variable behaviour depending 



FIG. 2: Largest Lyapunov exponent. A, is plotted as a func- 
tion of e for scale free network and f{x) = 4a;(l — x) and 
g{x) = X. Different regions are labeled as in Fig.0 



on e and initial conditions), the partially ordered region 
and the ordered region as shown by regions I to IV in 
Figs. (1) and (2). The symbols T, S, M, DQ, DP and 
F correspond to turbulent, self-organized, mixed, driven 
quasiperiodic, driven periodic and fixed point behaviour. 
For small values of e, we observe the turbulent behavior 
with all nodes evolving chaotically and there is no phase 
synchronization (region I-T). There is a critical value of 
coupling strength Cc beyond which synchronized clusters 
can be observed. This is a general property of all CMNs 
and the exact value of depends on the type of network, 
the type of coupling function and the parameter fj,. 

As e increases beyond Ec we get into a variable region 
(region II-S) which shows a variety of phase synchronized 
behavior, namely ordered chaotic, ordered quasiperiodic, 
ordered periodic and partially orderedbehaviour depend- 
ing on the initial conditions. The next region (region 
III-M) shows partially ordered chaotic behavior. Here, 
the number of clusters as well as the number of nodes 
in the clusters depend on the initial conditions and also 
they change with time. There are several isolated nodes 
not belonging to any cluster. Many of these nodes are of 
the floating type which keep on switching intermittently 
between an independent evolution and a phase synchro- 
nized evolution attached to some cluster. Last two re- 
gions (IV-DQ and IV-DP) are ordered quasiperiodic and 
ordered periodic regions showing driven synchronization. 
In these regions, the network always splits into two clus- 
ters. The two clusters are perfectly anti-phase synchro- 
nized with each other, i.e. when the nodes belonging to 
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FIG. 3: The figure shows several examples illustrating the self-organized and driven phase synchronization. The examples are 
chosen to demonstrate two diflerent ways of obtaining synchronized clusters and the variety of clusters that are formed. All the 
figures show node verses node diagram for N — Nc — 50. After an initial transient (about 2000 iterates) phaA'' — Nc = 50. Afte 
r an initial transient (about 2000 iterates) phase synchronized clusters are studied for T = 100. The logistic map parameter 
fi = 4 and coupstudied for T = 100. The logistic map parameter ^ = 4 and coupling function g{x) = x. The solid circles show 
that the two corresponding nodes are coupled and the open circles show that the corresponding nodes are phase synchronized. 
In each case the node numbers are reorganized so that nodes belonging to the same cluster are numbered consecutively and 
the clusters get displayed in decreasing sizes, (a) Figure shows turbulent phase for e = 0.10. (b) An ideal self-organized phase 
synchronization for e = 0.16. (c) Mixed behavior for e = 0.32. (d) A ideal driven phase synchronization for e = 0.90. The scale 
free networks were generated with No — 1 and m = 1. 



one cluster show minima those belonging to the other 
cluster show maxima. 

We now investigate the nature of phase ordering in 
different regions of the phase diagram. Fig. |31 shows 
node-node plots of the synchronized clusters with any 
two nodes belonging to the same cluster shown as open 



circles and the couplings between the nodes (Cy = 1) 
shown as solid circles. For small coupling strength, i.e. 
region I-T, nodes show turbulent behaviour and no clus- 
ter is formed (Fig. |3fa)). In region II-S the dominant 
behaviour is of self-organized type. Fig. |2Ib) shows an 
ideal self-organized synchronization with two clusters ob- 



6 



served in the middle of region II- S. Here, we observe that 
all the couplings except one are of intra-cluster type. Ex- 
actly opposite behavior is observed for the regions IV-DQ 
and IV-DP. Fig. Old) shows an ideal driven synchroniza- 
tion obtained in the middle of region IV-DP. Here, we 
find that all the couplings are of inter-cluster type with 
no intra-cluster couplings. This is clearly the phenomena 
of driven synchronization where the nodes of one cluster 
are driven into a phase synchronized state due to the 
couplings with the nodes of the other cluster. The phe- 
nomena of driven synchronization in this region is a very 
robust one in the sense that it is obtained for almost all 
initial conditions, the transient time is very small, the 
nodes belonging to the two clusters are uniquely deter- 
mined and we get a stable solution. In region III-M we 
get clusters of mixed type (Fig. OIc)), here the inter- 
cluster connections and the intra-cluster connections are 
almost equal in numbers. 

Mechanism of cluster formation: We observe that for 
small values of e the self organized behavior dominates 
while for large e driven behavior dominates. As the cou- 
pling parameter e increases from zero and we enter the 
region II-S, we observe phase synchronized clusters of the 
self organized type. Region III-M acts as a crossover re- 
gion from the self-organized to the driven behavior. Here, 
the clusters are of the mixed type. The number of inter- 
cluster couplings is approximately same as the number 
of intra-cluster couplings. In this region there is a com- 
petition between the self-organized and driven behavior. 
This appears to be the reason for the formation of several 
clusters and floating nodes as well as the sensitivity of 
these to the initial conditions. As e increases, we get into 
region IV-DQ where the driven synchronization domi- 
nates and most of the connections between the nodes are 
of the inter-cluster type with few intra-cluster connec- 
tions. This driven synchronization is further stabilized 
in region IV-DP with two perfectly anti-phase synchro- 
nized driven clusters. 

Quantitative measure for self-organized and driven be- 
haviour. To get a clear picture of self-organized and 
driven behaviour we define two quantities /intra and /inter 
as measures for the intra-cluster couplings and the inter- 
cluster couplings as follows: 



/intra 
/inter 



A^i, 



Ainter 



(6) 
(7) 



where Mntra and A^inter are the numbers of intra- and 
inter-cluster couplings respectively. In A^inter, couplings 
between two isolated nodes are not included. 

Fig.01shows the plot of /intra and /inter as a function of 
the coupling strength e. The figure clearly shows that for 
small coupling strength (region I-T) both /intra and /inter 
are zero indicating that there is no cluster formation at 
all, this is the turbulent region. As the coupling strength 
increases (e greater than some critical value ec) we get 
/intra ~ 1 at € ~ 0.2 (rcgion II-S). It shows that there 



are only intra-cluster couplings leading to self-organized 
clusters. As coupling strength increases further /intra de- 
creases and /inter increases i.e. there is a crossover from 
self-organized to driven behavior (regions III-M) . As cou- 
pling strength enters regions IV-DQ and IV-DP, we find 
that /inter is large which shows that in this region most 
of the connections are of the inter-cluster type. In region 
IV-DP we get /inter almost one corresponding to an ideal 
driven synchronized behaviour. 




0,4 - 



FIG. 4: The fraction of intra-cluster and inter-cluster cou- 
plings, /inter (solid clrcles) and fintra (open circle) are shown 
as a function of the coupling strength e for the scale-free net- 
works with f{x) — 4a;(l — x) and g{x) — x. The figure is 
obtained by averaging over 20 realizations of the network and 
50 random initial conditions for each realization. 

Behaviour of individual nodes forming clusters: Figs.El 
(a) and (b) show plot of time evolution of some typi- 
cal nodes. Fig. \^a) is for nodes in self-organized region 
(e = 0.15), where nodes belonging to the same cluster 
are marked with the same symbols. It is clearly seen 
that nodes with the same symbols i.e. belonging to the 
same cluster are phase synchronized and those belonging 
to different clusters are completely anti-phase synchro- 
nized, i.e. when the nodes in one cluster are showing 
minima, the nodes in other cluster are showing maxima. 
(This behaviour is observed for driven behaviour where 
two clusters are formed, i.e. nodes belonging to differ- 
ent clusters are anti-phase synchronized with each other.) 
Fig. I3b) plots the time evolution of three nodes in the 
partially ordered region (e = 0.35). We see that these 
nodes are not phase synchronized with each other. 

Now we explore different regions further to understand 
time evolution of individual nodes attached to some spe- 
cific cluster. Fig. plots all the nodes belonging to a 
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FIG. 5: Figures show time evolution of nodes belonging to different clusters. This figure is plotted for scale-free network with 
50 nodes coupled with g{x) = x. (a) A few nodes belonging to two phase synchronized clusters are shown. Nodes denoted by 
circles belong to one cluster and nodes denoted by squares to another cluster. Here e — 0.15. (b) Time series of three nodes 
which are not phase synchronized with each other arethree nodes which are not phase synchronized with each other are shown 
with three different symbols. Here, e — 0.35. 



cluster as a function of time, symbols indicate the time 
for which a given node belongs to a cluster. First we con- 
sider region II-S. Fig.EIa) shows a set of nodes (crosses) 
belonging to a cluster in the mixed region for e = 0.19 
and another set of nodes (open circle) belonging to an- 
other cluster for the same e but obtained with different 
initial conditions. For this e value all the nodes form self- 
organized clusters with no isolated nodes and this clus- 
ter formation is stable. It is seen from Fig. IHl^a) that all 
the nodes are permanent members of the cluster. Also 
comparing the members of two clusters which are ob- 
tained from different initial conditions we see that there 
are some common nodes while some are different. The 
reason is that these are self-organized clusters and this 
organization is not unique (see subsection 3. A. 4). On the 
other hand, driven synchronization (region IV-DP) leads 
to a unique cluster formation and does not depend on 
the initial conditions. 

Next we look at e = 0.4 (region III-M) where we get 
several clusters with some isolated nodes. In Fig. Elb), 
nodes belonging to a cluster are plotted as a function of 
time. We observe that there are some nodes which are 
attached to this cluster, intermittently leave the cluster, 
evolve independently or get attached with some other 
cluster and after some time again come back to the same 
cluster. These nodes are of the floating type which keep 
on switching intermittently between an independent evo- 
lution and a phase synchronized evolution with some 



cluster. For example, node number 12 in Fig. (tseries- 
clus)(b), which forms phase synchronized cluster with 
other nodes, in between leaves the cluster and evolve 
independently for some time. Time it spends with the 
cluster is about 90%. On the other hand node number 
24 evolves independently for almost 90% of the time and 
evolves in phase synchronization with the cluster for the 
rest of the time. 

Let T denote the residence time of a floating node in 
a cluster (i.e. the continuous time interval that the node 
is in a cluster). Fig. [7| plots the frequency of residence 
time /(r) of a floating node as a function of the residence 
time r. A good straight line fit on log-linear plot shows 
exponential dependence, /(t) ~ exp(— r/r^) where is 
the typical residence time for a given node. 



3. Nonlinear coupling 

Now we discuss the results for the nonlinear coupling 
of the type g{x) = f{x). Phase space diagram in the 
/i — e plane is plotted in Fig. |S1 Here we do not get 
clear and distinct regions as we get for g[x) = x form of 
coupling. Again the phase diagram is divided into dif- 
ferent regions I to V, based on the number of clusters as 
given in the beginning of this section. For /i < 3.5, we 
get coherent behaviour (regions V-P and VTF). To de- 
scribe the remaining phase diagram first consider /i = 4 
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FIG. 6: Figure shows the time evolution of nodes in a cluster for scale-free network, (a) shows two stationary clusters of 
self-organized type for e = 0.19 and g{x) = x. The two clusters are for the same e value but for two different initial conditions. 
The nodes belonging to the two clusters are denoted by open circles and crosses. Note that some nodes are common to both 
the clusters while some are different .nodes are common to both the clusters while some are different . This illustrates the 
nonuniqueness of nodes belonging self-organized clusters depending on the initial conditions, (b) shows a cluster with some 
permanent nodes and some floating nodes. Here e = 0.4 and g{x) — x. Node number 12, 24, 38 and 50 are oe = 0.4 and g{x) = x. 
Node number 12, 24, 38 and 50 are of floating type. They spend some time intermittently in a synchronized evolution with the 
given cluster and the remaining time in either a synchronized evolution with other clusters or in an independent evolution as 
an isolated node. 



line. Figure IHl shows the largest Lyapunov exponent as 
a function of the coupling strength e for /i = 4. For 
small coupling strengths no cluster is formed and we get 
the turbulent region (I-T). As the coupling strength in- 
creases we get into the variable region (II-D). In this re- 
gion we get partially ordered and ordered chaotic phase 
depending on the initial conditions. In a small portion 
in the middle of region II-D, all nodes form two ideal 
driven clusters. These two clusters are perfectly anti- 
phase synchronized with each other. Interestingly the 
dynamics still remains chaotic. In region III-T, we get 
almost turbulent behaviour with very few nodes form- 
ing synchronized clusters. Regions III-M and III-D are 
partially ordered chaotic regions. In these regions some 
nodes form clusters and several nodes are isolated or of 
the floating type. 

We now investigate the nature of phase ordering in dif- 
ferent regions of the phase diagram. Figs. ^1 are node- 
node plots showing different clusters and couplings (as 
in Fig. 13) for different e values belonging to different re- 
gions. In Fig. [rUr a) (region II-D) we observe that all the 
couplings are of inter-cluster type with no intra-cluster 
coupling. This is the phenomenon of driven synchroniza- 
tion. There are two clusters which are perfectly anti- 



phase synchronized. Fig. llUf b) is plotted for e in re- 
gion III-M and it shows clusters of different types. The 
fifth cluster has only inter-cluster couplings (driven type) 
while the remaining clusters have dominant intra-cluster 
couplings (self-organized type). There are several iso- 
lated nodes also. Fig. Iiur c) shows clusters where the 
driven behaviour dominates. It is interesting to note 
that for scale-free network and for this type of nonlinear 
coupling largest Lyapunov exponent is always positive 
(Fig.O i.e. the whole system remains chaotic but we get 
phase-synchronized behavior. 

Fig. II II shows the plot of /intra and /inter as a function 
of the coupling strength e for /i = 4. For small coupling 
strength both quantities are zero showing the turbulent 
region, and as the coupling strength increases clusters are 
formed, /inter is one at e w 0.13 which shows that here all 
the nodes are forming clusters and the clusters are of the 
driven type having only inter-cluster connections. As the 
coupling strength increases further, /inter and /intra be- 
come almost zero (region III-T) and subsequently start 
increasing slowly (region III-M) but we see that /inter 
is always greater than /intra leading to dominant driven 
phase synchronized clusters. For e > 0.7, /intra starts de- 
creasing and for e > 0.8, the driven behaviour becomes 
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FIG. 7: The figure plots the frequency of residence time /(t) 
of a floating node in a cluster as a function of the residence 
time r. The data is for node no 12 in Fig. |SIb). A good 
straight line fit on log-linear plot shows exponential depen- 
dence. 
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FIG. 9: Largest Lyapunov exponent, A, is plotted as a func- 
tion of e for scale-free network and f{x) — 4x{l — x) and 
g{x) = f{x). Different regions are labeled as in Fig. |H| 
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FIG. 8: Phase diagram showing turbulent, phase synchro- 
nized and coherent regions in the two parameter space of 
/.i and e for scale free network for j(x) = pix{\ — x) and 
g(x) = j(x). The determination of region boundaries and 
their classification and symbols are as explained in Fig. 
Calculations are for A'' = 50, m = 1,T = 100. The inset 
shows the phase diagram for the entire range of parameter /.i 
i.e. from to 4. 



more prominent (region III-D and Fig. llOf c)). For the 
regions III-M and III-D, we get phase synchronized clus- 
ters but the size of clusters as well as the number of nodes 
forming clusters both are small. 



J^. Network geometry and Cluster formation 

Geometrically, the organization of the scale-free net- 
work into connections of both self-organized and driven 
types is always possible for m = 1. For m = 1, our 
growth algorithm generates a tree type structure. A tree 
can be broken into different clusters in two distinct ways. 

(a) A tree can be broken into two or more disjoint clusters 
with only intra-cluster couplings by breaking one or more 
connections. Clearly, this splitting is not unique. This 
behaviour is observed in region II-S of the phase-diagram 
in Fig. 1 and can be seen by comparing Fig. 3(b) of this 
paper (two ideal self-organized cluster of sizes 41 and 9) 
and Fig. l(a)ofRef. [S^ (two ideal self-organized cluster 
of sizes 36 and 14) which are plotted for the same scale 
free network and g{x) = x but for different e values. 

(b) A tree can also be divided into two clusters by putting 
connected nodes into different clusters. This division is 
unique and leads to two clusters with only inter-cluster 
couplings. This behaviour is observed in region IV-DP 
of the phase-diagram in Fig. 1 and can be seen by com- 
paring Fig. 3(d) of this paper and Fig. 1(b) of Ref. [s^l 
which are again plotted for the same scale free network 
but for diferrent e values. 
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FIG. 10: The figure shows several examples illustrating the 
phase synchronization for scale-free network with coupling 
form g{x) = f{x) using node verses node diagram for — 
Nc = 50. After an initial transient (about 2000 iterates) 
phase synchronized clusters are studied for T — 100. The 
logistic map parameter /x = 4. The solid circles show that 
the two corresponding nodes are coupled and the open circles 
show that the corresponding nodes are phase synchronized. 
In each case the node numbers are reorganized so that nodes 
belonging to the same cluster are numbered consecutively and 
the clusters get displayed in decreasing sizes, (a) Figure show 
an ideal driven phase synchronization for e = 0.13. (b) Mixed 
behavior for e = 0.71. (c) A dominant driven behavior for 
e — 0.88. The scale free networks were generated with A'^o = 1 
and m = 1. 
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FIG. 11: The fraction of intra-cluster and inter-cluster cou- 
plings, finter (solid circles) and fintra (open circle) are shown 
as a function of the coupling strength e for the scale-free net- 
works with g{x) = f{x). The figure is obtained by averaging 
over 20 realizations of the network and 50 random initial con- 
ditions for each realization. 



For m > 1 and g(x) = x the dynamics of Eq. ^ leads 
to a similar phase diagram as in Fig. (1) with region 
II-S dominated by self-organized synchronization and re- 
gions IV-DQ and IV-DP dominated by driven synchro- 
nization. Though perfect inter- and intra-cluster cou- 
plings between the nodes as displayed in Figs. (3b) and 
(3d) are no longer observed, clustering in the region II-S 
is such that most of the couplings are of the intra-cluster 
type while for the regions IV-DQ and IV-DP they are of 
the inter-cluster type. As m increases the regions I and 
II are mostly unaffected, but the region IV shrinks and 
the region III grows in size. Fig. 1 121 shows different types 
of clusters in node-node diagram for different coupling 
strength. Fig. I12f a) is plotted for to = 3 in the variable 
region (e = 0.19) with nodes forming two clusters. It is 
clear that synchronization of nodes is mainly because of 
intra-cluster connections but there are a few inter-cluster 
connections also. Fig. ll2r bi is plotted for the ordered pe- 
riodic region at coupling strength e — 0.78, here the clus- 
ters are mainly of the driven type but they have intra- 
cluster connections also. In Figures I12f a'l and (b) the 
average degree of a node is 6, and breaking the network 
into clusters with only inter-cluster or intra-cluster cou- 
plings is not possible. As the average degree of a node in- 
creases further self-organized behaviour starts dominat- 
ing and for most of the e values nodes behave in phase- 
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synchronized manner forming one big cluster. 

For m > 1 and g{x) = f{x) we get similar kind of 
behaviour as for m — I with dominant driven clusters 
for most of the coupling strength region, but we do not 
get any ideal driven clusters. Fig. I12r c') is plotted for 
coupling strength e = 0.9 and m = 3. As m increases the 
region I showing turbulent behaviour remains unaffected, 
but the mixed region II grows in size while the region III 
shrinks. As m increases more and more nodes participate 
in cluster formation. The driven behaviour decreases in 
strength with increasing m and self-organized behaviour 
increases in strength. For m = 10, all nodes form one 
cluster for larger e values which is obviously of the self- 
organized type CFig. lTW dl). 

We have also studied the effect of size of the network 
on the synchronized cluster formation. The phenom- 
ena of self-organized and driven behavior persists for the 
largest size network that we have studied {N = 1000). 
The region II showing self-organized or driven behavior 
is mostly unaffected while the ordered regions showing 
driven behavior for large coupling strengths show a small 
shrinking in size. 



B. Coupled maps on one dimensional network 

For one dimensional CMN, each node is connected with 
m nearest neighbors (degree per node is 2m). First we 
consider m = 1, i.e each map is connected with just next 
neighbors on both sides. Fig ICT a) and (b) show /intra 
and /inter verses e for g{x) = x and g(x) — f{x) respec- 
tively and II = A, N ^ 50. For g{x) = x, after an initial 
turbulent region (e > ec), nodes form self-organized clus- 
ters (region II-S in Fig.^l and as the coupling strength in- 
creases we observe a crossover to driven clusters. The be- 
haviour of clusters as well as Lyapunov exponent graphs 
are similar to the scale-free network with the coupling 
form f{x) = x. Note that the nearest neighbor CMN 
with m = 1 is a tree and can be geometrically organized 
into both self-organized and driven type of clusters. 

For g{x) — f{x) coupling and m = 1, there is cluster 
formation for only small coupling strength region (corre- 
sponding to region II-D of Fig.lSJ as seen from Fig^Jb). 
Fig.ll4lshows largest Lyapunov exponent as a function of 
e for g{x) = f{x) and /i = 4.0. In region II-D the largest 
Lyapunov exponent is positive or negative, depending on 
the initial conditions and e values and for the rest of the 
coupling strength region Lyapunov exponent is positive. 
In region II-D synchronized clusters of driven type are 
seen. For larger coupling strengths where the system 
remains in chaotic zone, clusters are rarely formed and 
when formed are small in size. 

Fig. 1151 shows node-node plot showing synchronized 
clusters. Fig. I15f a) shows one self-organized cluster in 
region II-S (see Fig. for g{x) = x. In this region we 
also get two self-organized clusters depending on the ini- 
tial values and e. Fig. I15r b) shows two clusters of mixed 
type as well as several isolated nodes for e in region III- 



M for g{x) = x. Figures I15r c'l and I15r d') show driven 
clusters for g{x) — f{x) for e values in regions II-D of 
Fig.m 

We now consider the case to > 1. For g{x) — x we ob- 
serve self-organized clusters with some inter-cluster con- 
nections for the coupling strength region II-S and as the 
coupling strength increases there is a crossover to driven 
clusters. As the coupling strength increases further for 
e > 0.7 instead of forming driven clusters (as is observed 
for TO = 1) nodes form one synchronized cluster. As 
TO increases f intra increases and for e > we observe 
dominance of self-organized behaviour. For to = 5, and 
for coupling strength e > ec(« 0.13), all nodes form one 
or two clusters. For one cluster fintra = 1 and for two 
clusters intra-cluster and inter-cluster couplings are al- 
most equally distributed. For very large value of cou- 
pling strength (e > 0.7) we get clusters of dominant self- 
organized type. As the number of connections increases 
and typically becomes of the order that is a globally 
coupled state, we get one cluster of self-organized type. 

For g(x) = f{x) and m > 1, we find that as the num- 
ber of connections increases for small coupling strength 
(region II-D) we get two dominant driven phase synchro- 
nized clusters. For large coupling strength the number of 
nodes forming clusters and the sizes of clusters both in- 
crease with the increase in number of connections in the 
network. This behaviour is seen in Figs. I13r c) and (d) 
which show /intra and /inter verses e for g{x) = f{x), 
fi ^ A and respectively for m = 5 and to = 10. 
Fig. llGr a) shows the fraction of nodes forming clusters 
as a function of the number of connections Nc normal- 
ized with respect to the maximum number of connections 
Nm = N{N — l)/2 for two values of e. The overall in- 
crease in the number of nodes forming clusters is clearly 
seen. Fig. llBf b) shows the fraction of nodes in the largest 
cluster as a function Nc for two values of e. The overall 
growth in the size of the clusters with Nc is evident. 

Cluster formation with large number of connections (of 
the order of N"^) and its dependence on coupling strength 
is discussed in Refs. |^ [g^. It is reported that for 
these networks it is the coupling strength which affects 
the synchronized clusters and not the number of connec- 
tions. We find that when the number of connections is of 
the order of N there are significant deviations from this 
reported behaviour. We find that the size of the clus- 
ters and number of nodes forming clusters increases as 
the number of connections increase as discussed above. 
This behaviour approaches the reported behaviour as the 
number of connections increases and becomes of the or- 
der of N^. 



C. Coupled maps on Small world network 

Small world networks are constructed using the follow- 
ing algorithm by Watts and Strogatz 3]. Starting with 
a one-dimension ring lattice of N nodes in which every 
node is connected to its nearest k neighbors {k/2 = m 
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FIG. 12: Figure illustrates the cluster formation for the scale-free network as node vs node plot for A'" = 50 as in Fig. but 
with larger number of connections, (a) and (b) are plotted for g{x) — x and m = 3 and respectively for e = 0.19 and e = 0.78. 
(c) and (d) are plotted for g{x) — f{x), e — 0.90 and respectively for m = 3 and m — 10. 



on either side), we randomly rewire each connection of 
the lattice with probabiUty p such that self-loops and 
multiple connections are excluded. Thus, p = gives a 
regular network and p = 1 gives a random network. Here 
we present results for TV = 50 and m = 1. Figs. llTr a) 
andinib) plot /i„tra and /inter for g{x) = x and g{x) = 
f{x) respectively as a function of e for /i = 4. We find 
that for g{x) = x, the behaviour is very similar to that 



for the scale free networks and one-d lattice. We get self- 
organized clusters for e > Ec and there is a crossover to 
driven behavior as epsilon increases (Fis MYi a,) ). But for 
g{x) = f{x), nodes form clusters only for region 11-D of 
coupling strength and there is almost no cluster forma- 
tion for larger values of e ( Fig. llTf b) ). This behaviour 
changes as k increases and we observe some clusters for 
large e values also. Figure [THT a) shows node-node plot 



13 




FIG. 13: The fraction of intra-cluster and inter-cluster couplings, finter (solid circles) and fintra (open circles) are shown as a 
function of the coupling strength e. Figures (a) and (b) are for the one-d coupled maps with nearest neighbor coupling (m = 1) 
and for g{x) = x and g{x) = f{x) respectively. Figures (c) and (d) are for g{x) = f{x) and respectively m = 5 and m — 10. 
The figures are for A'^ = 50 and are obtained by averaging over 50 random initial conditions. 



of clusters for e = 0.45, m = 1 and g{x) = x showing 
dominant driven clusters. 



D. Coupled maps on Caley Tree 

We eenerate a Caley tree using the algorithm given in 
Ref. [23 ■ Starting with three branches at the first level, 
we split each branch into two at subsequent levels. For 
g{x) = X, the behaviour is similar to all other networks 
with the same number of connections (Fig. llTf c)). Fig- 
ure llSf b) shows node-node plot of two ideal driven phase 
synchronized clusters for e = 0.92, k = 2, N — 47 and 
g{x) = X. For g(x) = f{x) all nodes form driven clusters 
for region II-D, and for larger coupling strengths about 
40% of nodes form clusters of driven types (Fig. [T7jl (d)). 

E. Coupled maps on higher dimensional lattices 

Coupled maps on higher dimensional lattices also form 
synchronized clusters. First we give the result for two-d 



square lattices. Figs. lTW al and lTW b) plot /intra and /inter 
for g{x) — X and g{x) — f{x) respectively as a function 
of e for /i = 4. For g{x) — x the cluster formation is 
similar to other networks described earlier except for very 
large e close to one where we get a single self-organized 
cluster. For g{x) — f{x) cluster formation is similar to 
that in one-d networks with nearest and next nearest 
neighbor couplings. In small coupling strength region 
II-D (see Figure IHl , nodes form two clusters of driven 
type and for large coupling strength also driven clusters 
are observed with 25-30% nodes showing synchronized 
behaviour (Figure. I19f b)). Figures ITST c) and (d) show 
node-node plot of self-organized behaviour for g{x) — x 
and dominant driven behaviour for g{x) = f{x). 

Coupled maps on three-d cubic lattice (degree per node 
is six) for g(x) — x show clusters similar to the other 
networks discussed earlier. For g(x) — /(x), nodes form 
driven type of clusters at small coupling strength (region 
II-D) and mainly we observe three clusters. For large 
coupling strengths also nodes form driven clusters and 
the nodes participating in cluster formation is now much 
larger than the two-d case. 
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FIG. 14: Largest Lyapunov exponent, A, is plotted as a 
function of e for 1-d nearest neighbor coupled network for 
f{x) = 4a;(l — a;) and g{x) = f{x). Different regions are 
labeled as for scale-free network (see Fig. [Sjand Fig.|51) 



F. Coupled maps on random network 

Random networks are constructed by connecting each 
pair of nodes with probability p. First consider the case 
where the average degree per node is two. For Unear cou- 
pling g{x) = X cluster formation is the same as for other 
networks with same average degree. For g{x) = f{x) 
driven type clusters are observed in region II-D and no 
significant cluster formation is observed for larger cou- 
pling strengths. This behaviour is similar to one-d net- 
work with k = 2 but different from the corresponding 
scale free network. For coupled maps on random net- 
works with average degree per node equal to four and 
g{x) = f{x), clusters with dominant driven behaviour 
are observed for all e > Cc- 



G. Examples of self-organized and driven behavior 

There are several examples of self-organized and driven 
behaviour in naturally occurring systems. An impor- 
tant example in physics that includes both self-organized 
and driven behavior, is the nearest neighbor Ising model 
treated using Kawasaki dynamics. As the strength of 
the Ising interaction between spins changes sign from 
positive to negative there is a change of phase from 
a ferromagnetic (self-organized) to an antiferromagnetic 
(driven) behavior. In the antiferromagnetic state, i.e. 
driven behavior, the lattice spits into two sub-lattices 



with only inter-cluster interactions and no intra-cluster 
interactions. 

Several other examples are discussed in Ref. . 



IV. UNIVERSAL FEATURES 

We have discussed the behavior of evolution of coupled 
dynamical elements on different networks. For small cou- 
pling strengths unto a critical value Cc nodes show turbu- 
lent behaviour and for e > ec they form interesting phase- 
synchronized clusters. The critical value Cc depends on 
the type of network, the type of coupling function and 
parameter fj,. We find that clusters arc formed because 
of intra-cluster connections and/or intcr-clustcr connec- 
tions depending on the network, the type of coupling and 
the coupling strength. 

The linear coupling of type g(x) = x shows both 
types of clusters (self-organized and driven) and show a 
behaviour which is universal for all the networks that 
we have studied. For networks, with number of con- 
nections of the order of N, initially for a small range 
of coupling strength (region II-S in Figure 1) nodes 
form self-organized clusters and as coupling strength in- 
creases there is a crossover and reorganization of nodes 
to driven clusters. This behaviour is observed for all the 
networks that we have studied i.e. scale-free network, 
random network, small world network, Caley tree net- 
work, 1-d nearest neighbor coupled network, 1-d next 
to next neighbor coupled network, 2-d network. As the 
number of connections, Nc, increases, the driven be- 
haviour for large coupling strengths is suppressed and 
self-organized behaviour starts dominating. As Nc in- 
creases further a clear identification of the two mecha- 
nisms, self-organized and driven, becomes more and more 
difficult. As Nc becomes of the order of iV^, for very 
large coupling strengths we observe one spanning cluster 
of self-organized type. 

For coupling function g{x) — f{x) the formation of 
clusters and nature of clusters both depend on the type of 
network. Initially for a small range of coupling strength 
values (region II-D in Fig. ^ nodes form driven clusters 
for all networks but as coupling strength increases clus- 
ter formation and size of clusters both depend on the 
type of network. For networks with average connectivity 
per node equal to two, we find that for large coupling 
strengths number of nodes forming clusters and size of 
clusters both are small. For large coupling strengths 
where largest Lyapunov exponent is positive, less than 
50% nodes form clusters for scale-free network and Caley 
tree. For other types of networks cluster formation is not 
significant. As the total number of connections increases 
from Nc N to Nc ^ N'^, the number of nodes forming 
clusters as well as the size of clusters increase. Again as 
in the case of g{x) = x for very large coupling strengths 
one large spanning cluster is observed. 

It is interesting to note that nodes can form two or 
more stationary clusters even though coupled dynamics 
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FIG. 15: The figure illustrates the cluster formation for one-d nearest neighbor network using node-node plot as in Fig.|3 (a) 
and (b) are for g{x) — x and e — 0.16, and e — 0.30 respectively, (c) and (d) are for g{x) = f{x) and e = 0.13 and e — 0.15 
respectively. 



is in the chaotic regime. Here, a stationary cluster means 
that constituents of the cluster are unique and once nodes 
form a stationary cluster they belong to that cluster for- 
ever, that is the structure of the cluster does not depend 
on time. For fi — A, we get stationary clusters when 
the nodes form two clusters, but they are not stationary 
when they form three or more clusters. Note that three 
or more stationary clusters can be formed for /i < 4. 

For = 4 if the largest Lyapunov exponent is negative, 
the variables show periodic behaviour with even period. 



For fjL < A the periodic behaviour can have both odd and 
even periods. 

We observe ideal behaviour of both the types that is 
all the nodes forming driven clusters {f inter = 1) or all 
the nodes forming self-organized clusters [f inter = !)■ bi 
most cases where ideal behaviour is observed, the largest 
Lyapunov exponent is negative or zero giving stable clus- 
ters. However, in some cases ideal behaviour is also ob- 
served in the chaotic region. 
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FIG. 16: Figure (a) shows the fraction of nodes forming clusters as a function of the fraction of couphngs Nc/Nm where 
Nm = N{N — l)/2. The figures are plotted for 1-d coupled maps with g{x) — fix) and for e = 0.49 (closed circles) and e = 0.7 
(open circles). The results are for A'' = 50 and are obtained by averaging over ICQ random initial conditions. Figure (b) shows 
the fraction of nodes in the largest cluster as a function of Nc/Nm for e = 0.49 (closed circles) and e = 0.7 (open circles). Other 
parameters are same as above. 



V. CIRCLE MAP 

We have studied cluster formation by considering circle 
map as defining the local dynamics, given by 

f{x) = X + uj + {k/27T)sin{2TTx), (mod 1) 

Due to the modulo condition, instead of using the vari- 
able xt, we use a function of xt such as sin(7ra;t) satisfy- 
ing periodic boundary conditions to decide the location 
of maxima and minima which are used to determine the 
phase synchronization of two nodes (Eq. (|SJ|). With cir- 
cle map also we observe formation of clusters with the 
time evolution starting from initial random conditions. 
Here we discuss the results with the parameters of the 
circle map in the chaotic region {oj — 0.44 and fc = 6). 
For linear coupling g(x) — x and scale- free networks with 
m = 1, for small coupling strength nodes evolve chaot- 
ically with no cluster formation. As coupling strength 
increases nodes form clusters for 0.21 < e < 0.25. In 
most of this region the nodes form two cluster and these 
clusters are mainly of the driven type except in the initial 
part, e « 0.21, where self-organized clusters can be ob- 
served. As the coupling strength increases nodes behave 
in a turbulent manner and after e > 0.60 nodes form clus- 
ters of dominant driven type. Here the number of nodes 
forming clusters and the sizes of clusters, both are small. 
For the one dimensional linearly coupled networks, for 
linear coupling the nodes form phase synchronized clus- 



ters for coupling strength region 0.21 < e < 0.25. The 
clusters are mainly of the driven type except in the initial 
part, e « 0.21, where they are of the self-organized type. 
For large coupling strength they do not show any cluster 
formation. 

For g{x) — f{x) we found very negligible cluster for- 
mation for the entire range of the coupling strength for 
both scale free and one-d network. However, as m in- 
creases the nodes form phase synchronized clusters for e 
larger than some critical Cc- 

For the circle map the normalization factor (1 — e) in 
the first term of Eq. is not necessary and the following 
modified model can also be considered. 

N 

xl^, = f{xl) + ^Y. C^39{xi\ (mod 1). (8) 

'^^ 3 = 1 

We now discuss the synchronized cluster formation for 
the same parameter values as above [lu = 0.44 and fc = 6) 
for this modified model. For linear coupling, clusters are 
formed only for 0.02 < e < 0.17 with dominant self- 
organizied behaviour for most of the range except near 
e « 0.17 where the behaviour is of doninant driven type. 
For the scale free networks (m = 1) we have ordered 
states while for the one-d networks we have partially or- 
dered states. For nonlinear coupling, the clusters are 
formed for 0.0 < e < 0.09. The scale-free networks show 
mostly mixed type clusters while one-d networks show 
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FIG. 17: Prction of intra-cluster and inter-cluster couplings, finter (solid line) and fintra (dashed line) are shown as a function 
of the coupling strength e. Figures (a) and (b) arc for the small world network for g{x) — x and g{x) = j{x) respectively and 
iV = 50. Figures (c) and (d) are for the Caley tree with g(x) = x and g{x) = f{x) respectively and N = 47. The figures are 
obtained by averaging over 50 random initial conditions. Small world networks are generated with m = 1 and p = 0.06 [3]. 
Caley trees are generated with coordination number three [27]. 



dominant self-organized clusters. There is no cluster for- 
mation for larger coupling strengths for both linear and 
nonlinear coupling. However, as for the logistic map, syn- 
chronized clusters are obseved for large e as the number 
of connections inceases. 



VI. CONCLUSION AND DISCUSSION 

We have studied the properties of coupled dynamical 
elements on different types of networks. Wc find that in 
the course of time evolution they show phase synchro- 
nized cluster formation. We have mainly studied net- 
works with small number of connections {Nc ~ N) be- 
cause a large number of natural systems fall under this 
category of small connections. More importantly, with 
small number of connections, it is easy to identify the re- 
lation between the dynamical evolution, the cluster for- 
mation and the geometry of networks. We have studied 
the mechanism of cluster formation as well as behaviour 
of individual nodes either forming clusters or evolving 



independently. We use the logistic map for the local dy- 
namics and the two types of couplings. Depending upon 

the type of coupling, the regions for cluster formation 
and the behaviour of phase synchronization vary. We 
have identified two mechanism of cluster formation, self- 
organized and driven phase synchronization. 

By considering the number of inter- and intra-cluster 
couplings we can identify phase synchronized clusters 
with dominant self-organized behavior (S), dominant 
driven behavior (D) and mixed behavior (M) where both 
mechanisms contribute. Wo have also observed ideal clus- 
ters of both self-organized and driven type. In most 
cases where ideal behaviour is observed, the largest Lya- 
punov exponent is negative or zero giving stable clusters 
with periodic evolution. However, in some cases ideal be- 
haviour is also observed in the chaotic region. In most of 
the cases when synchronized clusters arc formed there are 
some isolated nodes which do not belong to any cluster. 
More interestingly there are some floating nodes which 
show an intermittent behavior between an independent 
evolution and an evolution synchronized with some clus- 




FIG. 18: Figure illustrates the cluster formation for different networks as node vs node plot as in Fig|3 (a) is plotted for small 
world network with e = 0.45, A'^ = 50 and g{x) — x. (b) is plotted for Caley tree with e = 0.92, N = 47 and g{x) — x. (c) and 
(d) are plotted for 2-d lattice (A^ — 49) with e = 0.21, g{x) = x and e = 0.19, g{x) = f{x) respectively. 



ter. The time spent by a floating node in the synchro- 
nized cluster shows an exponential distribution. 

By defining different states of the dynamical system 
using the number and type of clusters, we consider the 
phase-diagram in the /i — e plane for the local dynamics 
governed by the logistic map. When the local dynam- 
ics is in the chaotic region, for small coupling strengths 
we observe turbulent behaviour. There is a critical value 
Ec above which phase synchronized clusters are observed. 
For g(x) = x type of coupling, self-organized clusters are 



formed when the strength of the coupling is small. As the 
coupling strength increases there is a crossover from the 
self-organized to the driven behavior which also involves 
reorganization of nodes into different clusters. This be- 
haviour is almost independent of the type of networks. 
For non linear coupling of type g{x) — /(x), for small 
coupling strength phase synchronized clusters of driven 
type are formed, but for large coupling strength number 
of nodes forming cluster as well as size of cluster both are 
very small and almost negligible for many network. Only 
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FIG. 19: The fraction of intra-cluster and inter-cluster cou- 
plings, /inter (solid line) and fintra (dashed line) are shown 
as a function of the coupling strength e for two-d lattice. Fig- 
ures (a) and (b) are for g{x) — x and g{x) — f{x) respectively. 
The figures are for A*' = 49 and are obtained by averaging over 
50 random initial conditions. 



Denoting the number of elements of a set by |.|, let, 

(1) a = nMj nA4- 

(2) b = \M nJVk] - a. 

(3) c = IJVj nA4| - a. 

(4) d= IM nA/;-| - a. 

(5) e ^ \J^i \ ~ b - d ~ a. 

(6) \J\fj\ -c^d~a. 

(7) g^\Nk\-b-c-a. 
We have 

riik = a + b 
rijk = a + c 
Tij-i = a + d 



Hi 

nk 



a + b + d + e 
a + c + d + f 
a + b + c + g 



Consider the combination 



where 



X = 



riik 



for scale-free networks and Caley tree show some cluster 
formation for large coupling strengths. 

As the number of connections increases, most of the 
clusters become of the mixed type where both the mech- 
anisms contribute. We find that in general, the self- 
organized behaviour is strengthened and also the number 
of nodes forming clusters as well as the size of clusters 
increase. As the number of connections become of the 
order of N"^ , self-organized behaviour with a single span- 
ning cluster is observed for e larger than some value. 

It is interesting to consider the dynamical origin of the 
self-organized and driven phase synchronization. A clue 
can be obtained by considering small networks of two 
or three nodes and also tailormade large networks such 
as globally coupled networks and complete bipartite net- 
works. These are considered in Ref. l66|. These studies 
reveal that the intra-cluster coupling term between the 
two varibles a;* and x-' , adds a decay term to the dynamics 

of the difference variable x'^ = x'^ — x^ leading to synchro- Case B. rii < Uk < n 
nization of these two variables. On the other hand, the 
inter-cluster coupling term between the two varibles a;* 

and x^ , cancels out in the dynamics of x'^ and the two X = 

variables now belong to different synchronized clusters. 

In this paper we have discussed the numerical results of < — 

phase synchronization on CMN and the two mechanisms 

of cluster formation. In another paper we discuss the < 1 

stabilty analysis of synchronized clusters in some simple 

networks that illustrate the two mechanisms of synchro- Case C. Uk < fii < Uj-. 
nization 1661. 



1-X 



(Al) 



rhk 



max(ni,nfc) max(nj,nfe) max(ni,nj) 



The triangle inequality is proved if A" < 1. Consider the 
following three general cases. 

Case A. rii < rij < Uk- 



X 



nk 



nk 



< 



< 1 



nk 



(A2) 



nk 



b + c 



nk 



(A3) 



APPENDIX A 



Here we show that the definition ||SJ) of phase distance 
dij between two nodes i and j satisfies metric properties. 
Let J\fi denote the set of minima of the variable in a 
time interval T. The phase distance satisfies the following 
metric properties. 
f^A.^ d^j — dji . 

(B) = only if M = A/'j- 

(C) Triangle inequality: Consider three nodes i, j and k. 



X = 



< 



< 



a + b ^ a + c a + d 

rii nj rij 

a + b + c 

a + b + c 



nk 



< 1 



(A4) 



This proves the triangle inequality. 
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